RL-Tfl-84-52 
In-HoiiM  Report 
<Ium1904 


AD-A281  188 


MONOPULSE  TRACKER  OPERATING 
CHARACTERISTICS  MEASUREMENT 
ERROR  PERFORMANCE  FOR  AZIMUTH 
TRACKING 


Scott  Bolen 


DTIC 

Selecte 

JUL  081994|  I 


APPROVED  FOR  PUBLC /PLEASE;  DISTfKBUTK)NUNLmTED. 


94-20710 

illllilli  . 


'^d 


QUALirr 


Aim 


Rome  Laboratory 
Air  Force  Materiel  Command 
Griffiss  Air  Force  Base,  New  York 


94 


154 


This  rsport  has  basn  rsvlswed  by  the  Rome  Laboratory  Public  Affairs  Office 
(PA)  and  Is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RL-TR-94-52  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


JOSEPH  J.  SIMONS,  Chief 


Wide  Area  Radar  Surveillance  Division 
Surveillance  &  Photonics  Directorate 


FOR  THE  COMMANDER: 


LUKE  L.  LUCAS,  Colonel,  USAF 
Deputy  Director 

Surveillance  &  Photonics  Directorate 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  Rome  Laboratory 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization 
please  notify  RL  ( OCDS  )  Griff iss  AFB  NY  13441.  This  will  assist  us  in  maintalnine 
a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


1 .  AOENCY  USE  ONLY  (Lmm  8lanl4 


2.  REPORT  0A1E 
Jua«  1994 


4.TTTlEAN0SUBTm£ 

MONOPULSE  TRACKER  OPERATING  CHARACTERISTICS  MEASUREMENT 
ERROR  PERFORMANCE  FOR  AZIMUTH  TRACKING 


&AUTN0R(9 

Scocc  Bolen 


a  REPORT  TVPE  AND  OATES  COVERED 
In-House  May  93  -  Jul  93 


a  PUNOMQ  NUMBERS 

PE  -  62702F 
PR  -  4506 
TA  -  14 
WU  -  TK 


7.  PERFORMMQ  OROAMZATION  NAME(S)  AND  AOORESS<ES) 
Rome  Laboratory  (OCOS) 

26  Electronic  PIcy 
Griffiss  AFB  NY  13441-4514 


a  PERFORMMQ  ORGANIZATION 
REPORT  NUMBER 
RL-TR-94-52 


a  SPONSOHMQMONrrOHMQ  AGENCY  NAiME<S)  AND  ADDRESSEES) 
Rome  Laboratory  (OCDS) 

26  Electronic  Pky 
Griffiss  AFB  NY  13441-4514 


ia  SPONSORMOMONTTORMC 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

Rome  Laboratory  Project  Engineer:  Scott  Bolen/OCOS  (315)  330-4441 


12a.  OtSTHBUTION/AVAtABUTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


lisa  DSmBURON  CODE 


1  a  ABSTRACTpmmrun  mD-aee 

This  paper  derives  the  operating  characteristics  for  a  monopulse  radar  tracker  as  a 
function  of  probability  of  detection  and  probability  of  false  alarm.  Measurement 
error  performance  is  developed  for  azimuth  tracking  in  an  amplitude  comparison  mono¬ 
pulse  system.  A  set  of  receiver  operating  characteristics  is  also  presented  to 
demonstrate  how  to  determine  the  detection  threshold  in  the  signal  processing  system 
to  minimize  tracking  errors  for  a  given  signal- to-noise  ratio. 


14.  SUBJECT  TERMS 

Monopulse  Tracker,  Tracker  Operating  Characteristics 


IS  numnen  of  poces 
20 _ 

tepncecooe 


1 7.  SECURTTY  CtASST  ICARON 
®*'l^^SIFIED 


1&  SECURTTY  CLASSMCATION  IQ  SECURTTY  CLASSFKATION  2a  UMTTATION  OF  ABSTRACT 
^^iM^FIED  ^W^UIEU  U/L 


Monopulse  Tracker  Operating  Characteristics 
Measurement  Error  Performance  for  Azimuth  Tracking 


Scott  Bolen 
July  1993 


Abstract:  This  paper  derives  the  operating  characteristics 
for  a  monopulse  radar  tracker  as  a  function  of  probability 
of  detection  and  probability  of  false  alarm.  Measurement 
error  performance  is  developed  for  azimuth  tracking  in  an 
amplitude • comparison  monopulse  system.  A  set  of  receiver 
operating  characteristics  is  also  presented  to  demonstrate 
how  to  determine  the  detection  threshold  in  the  signal 
processing  system  to  minimize  tracking  errors  for  a  given 
signal-to-noise  ratio. 


I .  Introduction 

The  amplitude  comparison  monopulse  system  employs  two 
overlapping  antenna  beams  radiating  simultaneously  to  obtain 
an  angular  error  in  one  coordinate.  The  monopulse  tracker 
is  essentially  a  closed-loop  control  system  that  follows  the 
change  in  position  of  a  moving  target.  Target  position  is 
indicated  by  the  relative  amplitude  of  the  receive  signals 
from  each  of  the  two  beams.  Figure  1  shows  the  monopulse 
measurement  process.  In  figure  1(a)  the  relative  beam 
pattern  is  shown.  There  are  two  symmetric  beams  that  can  be 
considered  to  be  offset  by  a  distance  of  +/-  x  along  an 
arbitrary  axis  X.  The  response  of  each  beam  can  be  defined 
as  a  function  f(y)  such  that: 

fl(y)  =  the  response  of  the  received  signal 

from  beam  one 

and 

f2(y)  =  the  response  of  the  received  signal 

from  beam  two. 

In  the  monopulse  system  the  receive  signal  from  each  beam  is 
subtracted  to  form  a  difference  response  or  error  signal  as 
shown  in  figure  1 (b) .  This  error  signal  is  given  by: 

error (y)  =  fl (y)  -  f2{y) .  (1) 
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The  difference  response  forms  the  control  signal  for  the 
closed-loop  monopulse  tracking  system  [1/2]. 

In  the  monopulse  system  a  null  is  formed  in  the  antenna 
pattern  between  the  two  radiated  beams.  The  focus  of  this 
paper  is  to  develop  a  set  of  operating  characteristics  for  a 
tracking  system  that  keeps  the  target  within  the  null  of  the 
antenna  pattern.  This  type  of  system  is  called  a  "null" 
tracker.  The  null  used  for  tracking  occurs  in  the  center  of 
the  antenna  or  array.  The  amplitude  of  the  error  signal 
specifies  the  degree  to  which  the  antenna  is  centered  on  the 
target  [1/2]. 


II.  Development  of  the  Monopulse  Tracker  Operating 

Characteristics 

The  difference  response  represented  in  figure  1(a) 
shows  that  at  the  null  of  the  antenna  pattern 

error (y)  =  0. 

In  an  unambiguous  track  condition/  this  response  indicates 
that  the  target  is  situated  in  the  center  of  the  antenna 
between  the  two  transmitted  beams.  A  positive  value  for  the 
error  signal  indicates  that  the  target  is  in  beam  one  and  a 
negative  value  places  the  target  in  beam  two. 

The  magnitude  of  the  difference  response  is  dependent 
on  the  received  signal  strength  of  the  target  and  the 
location  of  the  target  with  respect  to  the  location  of  the 
transmitted  beams  (ie  the  tracking  null) .  In  practical 
radar  systems/  the  antenna  pattern  developed  by  the 
transmission  of  the  two  beams  will  contain  other  nulls  in 
addition  to  the  desired  tracking  null.  The  presence  of 
these  unwanted  nulls  can  cause  ambiguities  to  occur  in  the 
tracking  process.  That  iS/  when  error (y)  =  0  this  can 
indicate  that;  (a)  the  target  is  contained  within  the 
tracking  null  or  (b)  the  target  is  in  another  null  off  angle 
from  the  primary  tracking  null.  The  standard  approach  to 
overcome  this  ambiguity  is  to  from  another  response 

sum(y)  =  fl(y)  +  f2(y)  .  (2) 

This  function  is  referred  to  as  the  Siam  response.  Sum(y)  is 
used  for  target  detection  and  can  also  be  used  to  avoid 
tracking  ambiguities  in  azimuth  measurements.  A  normalized 
difference  response  can  be  formed  given  by  the  function 
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enorni{y) 


error (y) 
suni{y) 


(3) 


fl(y)  -  f2(y) 

“  fl(y)  +  f2(y)  • 

Ambiguities  can  be  avoided  if  the  sum(y)  is  set  above  a 
given  detection  threshold. 

A.  Measurement  Error  in  Noise 

In  the  difference  pattern  response  shown  in  figure 
1 (b) ,  the  slope  of  the  error  response  as  it  crosses  through 
the  zero  point  on  the  measurement  axis  is  called  the 
difference  slope  of  the  monopulse  measurement  [2] .  The  rate 
of  change  in  the  slope  of  the  curve  at  this  point  indicates 
the  relative  measurement  sensitivity  of  the  system.  A 
sharply  rising  slope  indicates  a  highly  sensitive  system  and 
a  slow  rising  slope  indicates  a  less  sensitive  system.  In 
this  case  sensitivity  refers  to  the  dynamic  response  of  the 
tracking  system.  Barton  [2]  defines  a  normalized  difference 
slope  as  a  differential  function  given  by 

_  d(enorm(y) ) 

^  d(bmw/bmw3) 

where 

bmw  -  the  antenna  beamwidth, 

bmw3  -  the  half  power  antenna  beamwidth. 

*  bmw/bmw3  is  a  normalized  beamwidth  parameter 

In  the  presence  of  noise,  spurious  error  signals  can  be 
generated  that  can  corrupt  the  tracking  process.  It  is 
assumed  that  the  noise  is  a  random  variable  with  an  even 
probability  density  function  (pdf) .  For  additive  noise  such 
that 


Total 

Received 

Signal 


Received  +  noise 

Target 

Signal 


the  mean  of  the  noise  pdf  can  be  considered  to  be  centered 
at  the  actual  target  position.  For  this  case,  the  azimuth 
measurement  position  error  will  vary  symmetrically  about  the 
actual  target  position  on  the  measurement  axis  with  an  rms 
error  given  by  [2]  i  Distibution/ 
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Dist 


Avail  and/or 
Special 


sigma 


_ bmw3 _ 

km[2  (SNR)n]l/2 


(4) 


where  SNR  =»  signal-to-noise  ratio 

n  *  number  of  pulses  integrated. 

Equation  (4)  is  the  fundamental  relationship  that  will  be 
used  throughout  the  paper  to  describe  the  rms  position  error 
of  the  monopulse  estimate  in  a  thermal  noise  environment. 

B.  Measurement  Error  Due  to  Target  Dynamics 

In  classical  control  system  theory,  the  lag  error 
exhibited  by  a  closed-loop  tracking  system,  such  as  the 
monopulse  tracker,  can  be  described  in  terms  of  the  dynamics 
of  the  system.  Servomechanism  control  theory  can  be  used  to 
predict  the  lag  error.  In  an  established  formulation,  the 
lag  error  of  the  system  is  given  by  [2,3] 

wa  d(wa)  (wa) 

®  *  Kv  Ka  K3  +  •  •  • 


where  Kv,  Ka,  and  K3  are  servo  error  coefficients. 


In  this  paper,  it  is  assumed  that  the  target  acceleration 
contributes  significantly  to  the  dynamics  of  system  and  that 
the  other  kinematic  properties  of  the  target  are  negligible. 
In  this  case,  the  system  servo  error  (dynamic  lag)  can  be 
defined  simply  as 


ea 


d(wa) 
Ka  * 


(6) 


Barton  [2]  describes  the  servo  error  as 
~  (2.5)RBn2 

where  a  =  target  acceleration 

R  =  range  to  target  from  the  radar 

Bn  =  servo  bandwidth  (eq  noise  bandwidth) . 

The  variance  of  the  azimuth  measurement  error  specified 
as  a  function  of  thermal  noise  and  dynamic  lag  can  be 
determined  by  combining  equation  (4)  and  equation  (7)  (after 
Barton  [2]) 

var  =  (sigma) 2  +  (ea)2  (8) 
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_ (bmw3^)  Bn _  _ a£ _ 

{km2) (fr) (Bt) (SNR)  *  6.3(R2){Bn2) 

where  fr  =  pulse  repetition  frequency 
Bt  =  time-bandwidth  product. 

C.  Measurement  Error  Expressed  in  Terms  of  Pd  &  Pf 

The  signal-to-noise  ratio  (SNR)  refers  to  the  power 
ratio  of  the  detected  signal  to  noise  at  the  input  of  the 
radar  receiver.  Another  parameter  called  the  detectability 
factor,  D,  can  be  defined  in  the  same  way  as  the  SNR  except 
that  we  will  let  D  be  the  required  value  of  the  power  ratio 
calculated  from  the  probability  of  detection  (Pd)  and 
probability  of  false  alarm  (Pf) .  For  a  Swerling  Case  1 
fluctuating  target/  the  single  pulse  detectability  factor  is 
specified  by  [2] 


Dl(l) 


In(Pf) 
In (Pd) 


Substituting  equation  (9)  into  equation  (8)  yields 


var  = 


_ (bmw3) 2 Bn _ 

(km)2(fr) (Bt) (ln(Pf)/ln(Pd)  -  1) 


6.3(R)2(Bn)2* 

Equation  (10)  now  defines  the  variance  of  the  azimuth 
traclcing  error  in  terms  of  Pd  and  Pf.  Simplifying  and 
solving  for  Pd  yields 


exp ( (var  - 


(a)2 

6.3(R)2(Bn)2' 

_ ()an)2(fr)  (Bt)  (In(Pf) ) 


Sbmw3)2(Bn)  +  (var)  ()aa)  2  ( fr)  (Bt)  +  (a)  2  (fr)  (Bt)  ^  ^  ’ 


From  equation  (11)  the  Monopulse  TracJcer  Operating 
Characteristics  (MTOC)  for  single  pulse  detection  and  for  a 
Swerling  Case  1  target  can  be  calculated.  Contours  of  the 
measurement  variance  error  can  be  calculated  over  the  Pd-Pf 
plane. 
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Example  -  Computation  of  MTOC  for  a  Given  Set  of  Radar 
Parameters  (see  figure  2) : 

For  the  parameters: 

km  =  1 . 6 

fr  =  1000.0  Hz 

Bt  -  1.0 

R  =  100.0  nmi 

Bn  =  3.0  Hz 

a  =20.0  M/sec 

bmw3  =  0.02  radians. 

Lines  of  constant  azimuth  tracking  variance  errors  can  be 
computed  as  a  function  of  Pd  and  Pf.  Figure  2  shows  the 
monopulse  azimuth  variance  errors  or  precision  of  the 
monopulse  estimate  expressed  in  decibels.  Each  line  shown 
in  the  figure  forms  a  contour  for  a  specific  error.  The 
entire  set  of  contours  forms  the  MTOC  for  the  example 
problem. 


III.  Calculation  of  Receiver  Operating  Characteristic 

The  Receiver  Operating  Characteristic  (ROC)  can  be 
calculated  in  terms  of  the  Pd  and  Pf.  There  are  several 
examples  in  the  literature  that  describe  the  detection 
process  and  demonstrate  how  to  determine  the  ROC  [4] .  In 
this  paper  it  is  assumed  that  the  ROC  is  determined  from  a 
binary  hypothesis  test  that  forms  a  likelihood  ratio  such  as 
the  Neyman-Pearson  decision  test. 

In  the  likelihood  ratio  test  there  are  two  hypotheses 
that  are  formed  in  the  decision  process.  The  observations 
for  the  two  hypotheses  are 

Ho:  si  =  ni  i  =  1,2,...,N 

(a  target  is  not  present) 

and 

HI:  si  =  m  +  ni  i  =  1/2, ...,N 

(a  target  is  present) 

where  s  is  the  total  signal  taken  at  the  receiver  which 
is  equal  to  the  radar  return  signal  plus  noise;  m  =  target 
signal/  ni  =  noise  signal,  N  =  number  of  samples  taken. 
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Given  the  above  hypothesis,  the  probability  density  of  si 
can  be  described  as  a  set  of  conditional  probabilities  for 
each  hypothesis  such  that 


P  si/Hl(Si/Hl)  *  pni(Si  -  m) 

and 

P  si/Ho (Si/Ho)  =  pni{Si). 

For  noise  samples  that  are  Gaussian  with  zero  mean  and 
variance  of  (sigma) 2,  the  probability  density  function 
of  si  under  each  hypothesis  can  be  described  as 


pni (Si  -  m) 


and 


pni (Si) 


1 

yj2  (pi)  (sigma) 


exp  (- 


(Si  -  m)2^ 
2  (sigma)  2 


1 

yj2  (pi)  (sigma) 


exp  (- 


2 


(Si)2  ^ 

(sigma) 2 


Because  the  ni  are  statistically  independent,  the  joint 
probability  density  of  each  si  is  the  product  of  the 
individual  probability  densities  such  that 


P  s/Hl(S/Hl) 


and 


1 

•v/2  (pi)  (sigma) 


exp  (- 


(Si  -  m)2^ 
2 (sigma) 2 


P  s/Ho(S/Ho) 


W9 

TT 


“^2  (pi)  (sigma) 


exp  (- 


(Si) 


2 (sigma) 2 


)  . 


The  lilcelihood  ratio  is  given  by 


1  .  (Si  -  m)2. 

~r===—==  exp ( -  - o) 

V2 (pi) (sigma)  2 (sigma) ^ 

1  ^  (Si)2 

>/2  (pi)  (sigma)  2  (sigma)  2 


So  for  a  given  threshold,  the  decision  criteria  is 


R(S)  > 


threshold 


choose  HI 


R(S)  < 


threshold 


choose  Ho. 


Figure  3  shows  the  density  functions  for  calculating 
the  probability  of  detection  and  false  alarm  given  a 
Gaussian  distribution.  The  decision  threshold  is  marked  on 
the  figure  and  is  used  to  calculate  the  values  of  Pd  and  Pf 
where  [4] 

_ .In (threshold)  d. 

Pf  =  1  -  erf( - 5 -  +  2) 

and 


_ .In (threshold)  d. 
1  -  erf  ( - 5 - j) 


where  erf  denotes  the  error  function  and  d  is  the 
distance  between  the  means  of  the  two  distributions. 

Figure  4  shows  the  ROC  as  a  function  of  Pd  and  Pf. 

Each  line  forms  a  contour  of  the  signal-to-noise  ratio 
present  at  the  receiver  for  a  given  Pd  and  Pf.  The  entire 
set  of  contours  forms  the  ROC  for  the  given  hypothesis  test. 


IV.  Summary  -  Characterizing  System  Performance 

System  performance  can  be  completely  described  by  the 
probability  of  detection  and  false  alarm.  From  the  example 
problems,  the  system  performance  can  be  characterized  by 
overlaying  the  ROC  contours  onto  the  MTOC  contours  as  shown 
in  figure  5.  The  intersection  of  the  contours  defines  the 
operating  point  of  the  system.  For  a  given  Pd  and  Pf,  the 
receiver  signal-to-noise  ratio  can  be  found  and  the 
subsequent  tracker  error  can  be  determined  from  the  graph. 
Figure  5  can  be  used  to  either  predict  track  performance  or 
to  specify  receiver  performance  based  on  a  required  track 
performance. 

Also,  since  Pd  and  Pf  are  a  function  of  the  detection 
threshold,  a  particular  receiver  operating  point  can  be 
determined  by  specifying  a  detection  threshold.  This 
threshold  can  be  determined  by  trading  off  performance  in 
terms  of  Pd  and  Pf  which  also  defines  the  track  error. 
Hence,  system  performance  can  be  completely  defined  by  the 
value  of  the  detection  threshold  level. 
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Figure  1 


(a) 
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Monopulse  measurement  process:  (a)  monopulse  response, 
(b)  difference  response. 
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Probability  of  Detection  (Pd) 


1 


Probability  of  False  Alarm  (Pfa) 

Figure  2:  Monopulse  Tracker  Operating  Characteristic  (MTOC). 

Variance  of  tracking  error  in  azimuth  due  to  dynamic  lag 
and  thermal  noise. 
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Figure  3 


(a) 


Pd 


(b) 


Error  probabilities,  Gaussian  distribution:  (a)  Pfa  calculation, 
(b)  Pd  calculation. 
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Probability  of  False  Alarm  (Pfa) 

Figure  5:  Monopulse  tracker  system  performance  for  a  given  receiver 
operating  point. 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  sdencs  and 
technologies  of  (xmimand,  control,  communicatiortt  and  intefligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  deveiopment  and  test  pograms  in  all 
applicable  technologies; 

b.  Trar»itions  technology  to  current  and  futisB  systems  to  improve 
operational  capability,  readiness.  arKi  supportabffity; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
CommarKl  product  centers  and  other  Ah’  Force  organizations; 

d.  Promotes  transfer  of  techndogy  to  the  private  sector; 

e.  MairTtara  leading  edge  technolo^cal  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  inteHigenca,  reliability 
sdenca,  electro-magnetic  technology,  photonics,  signal  procrasing.  and 
computational  science. 


The  thrust  areas  of  technical  competence  indude:  SurveHiance, 
Communications,  Command  and  Control,  Intelligence,  Signal  f^ocessing. 
Computer  Sdence  and  Technology,  Electromagnetic  Technology. 
Photonics  and  Reliabiiity  Sciences. 


